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EVAPORATION  RATES  OF  CHEMICAL  WARFARE  AGENTS 
USING  5  CM  WIND  TUNNELS 
IV.  VX  FROM  GLASS 


1.  INTRODUCTION 

The  evaporation  rates  of  the  chemical  warfare  agent  sulfur  mustard  from  glass, 
concrete,  and  sand  as  a  function  of  temperature,  drop  size,  and  air  flow  rate  have  been  previously 
reported.113  Select  evaporation  and  degradation  rates  of  ton  container  VX  on  sand  and  glass 
have  also  been  reported  in  the  literature.4  In  this  work,  additional  evaporation  data  for  VX  on  the 
same  glass,  as  used  previously,  as  a  function  of  temperature,  drop  size,  and  air  flow  rate 
measured  in  5  cm  wind  tunnels  are  presented.  This  report  describes  the  data  analysis  and 
demonstrates  the  robustness  of  the  set  of  data  that  will  be  passed  to  the  modelers  for  eventual 
incorporation  into  field  models  such  as  VLSTRACK,  which  predicts  vapor  concentration  as  a 
function  of  environmental  conditions. 


2.  EXPERIMENTAL  PROCEDURES 

2  1  Wind  Tunnel 

The  5  cm  wind  tunnels  that  were  used  in  the  present  investigations  ha\e  been 
previously  described  and  were  tbe  same  as  those  used  for  similar  earlier  studies  on  glass.3  The 
design  of  the  5  cm  laboratory-sized  wind  tunnels6  and  the  wind  tunnel  characteristics  compared 
to  other  w  ind  tunnels  and  outdoors  measurements  have  been  published.  In  order  to  expose  the 
agent  to  the  wind  flow,  the  piston  was  removed  and  the  test  substrate  (a  1.5  in.  diameter  circle) 
with  the  droplet  of  agent  on  it  was  placed  onto  the  piston  and  inserted  into  the  wind  tunnel.  The 
humidified,  temperature-controlled  air  from  a  Miller-Nelson  Environmental  Control  Unit  (ECU) 
(tunnel  a)  or  an  Aalborg  MFC  (tunnels  e,  d,  k,  1)  was  then  passed  over  the  sample  and  the  vapors 
were  collected  on  Markes  Tenax'  thermal  desorption  tubes  (Agilent  Technologies,  Santa  Clara, 
CA)  at  the  vapor  sampling  inlet.  The  amount  of  agent  on  each  tube  was  measured  based  upon  a 
standard  in  the  Gas  chromatography/Mass  Spectrometry  (GC/MS).  The  sample  volume  and 
tunnel  air  flow  rate  were  known;  thus,  the  agent  concentration  (mg/m3)  and  evaporation  rate 
(pg/tnin)  could  be  calculated.  The  rates  were  not  calculated  for  the  initial  5  min  of  the 
experiment  (before  the  instrumentation  reached  equilibrium)  nor  at  the  end  of  the  experiment 
(when  concentration  of  mustard  was  nearing  a  plateau  due  to  sample  exhaustion),  lienee,  the 
data  in  the  middle  of  eaeh  experimental  run  were  used  to  calculate  the  evaporation  rates. 

Air  flows  were  18,  181,  and  405  standard  liters  per  minute  (SLPM),  which 
corresponded  to  velocity  values  at  a  1  cm  height  of  0.22,  1.7,  and  3.6  m/s.  The  flow  volume  per 
thermal  desorption  tube  was  typically  2  to  10  L  volume,  and  tbe  tubes  were  automatically 
switched  using  a  proprietary  Versatile  Tube  Sampler.  The  rate  at  w  hich  the  tubes  w  ere  switched 
w  as  adjusted  based  upon  the  evaporation  rate  of  the  agent.  The  air  and  substrate  temperatures 
investigated  in  this  work  were  35,  42,  and  50  °C,  and  the  droplet  sizes  were  1, 6,  and  9  pL, 
corresponding  to  droplet  diameters  of  1 .24,  2.25,  and  2.58  mm  and  contamination  densities  of 
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approximately  1.3,  7,  and  1 1  g/m  .  The  droplets  masses  used  in  the  caleulations  and  tables  were 
based  upon  the  pipette  setting;  the  samples  were  not  weighed. 

In  general,  the  evaporation  of  the  VX  was  measured  until  no  further  decrease  in 
the  concentration  of  vapor  was  detected  at  which  time  the  experiment  was  terminated. 
Summation  of  the  concentrations  over  time  yielded  the  cumulative  amount  of  VX  that  had 
evaporated.  The  data  collected  arc  shown  in  the  Appendix. 

2.2  Detection  of  VX  in  the  Wind  Tunnel  Effluent 

For  the  analysis  of  VX,  a  silver  fluoride  pad  (CAMSCO,  Houston,  TX)  was 
inserted  onto  the  end  of  Tenax  TA  thermal  desorption  tubes  (Markcs  International,  Llantrisant, 
UK)  to  convert  any  VX  in  the  stream  to  its  G-analog,  ethyl  methylphosphonofluoridate  (EMPF, 
EA-1 207)  which  is  more  volatile  and  therefore  easier  to  analyze  by  thermal  desorption  GC.  The 
chemical  equation  for  the  reaction  is  shown  in  Scheme  1 . 


O  C3H7 

H3C— P — S— CH2CH2N  +  AgF 
OC2H5  c3h7 


O  C3H7 

->  H3C — P — F  +  AgSCH2CH2N 


oc2h5 


I 

c3H7 


Scheme  1 .  V-to-G  Conversion  on  Silver  Fluoride  Pads. 


The  ethyl  methylphosphonofluoridate  was  desorbed  from  the  thermal  desorption 
tubes  using  a  Markes  UNITY/ULTRA  Desorption  system  (Markcs  International,  Llantrisant, 
UK),  and  analyzed  on  an  Agilent  6890/5973  Gas  Chromatography/Mass  Spectrometry  Detector 
(GC/MSD  [Agilent  Technologies,  Santa  Clara,  CA])  using  a  HP-5MS  capillary  column  (30  m 
long,  0.25  mm  id,  0.25  pm  film  thickness,  (5%-phenyl)-methylpolysiloxane  stationary  phase).9 
Oven  parameters  were  60  ( 1 .5  min)  to  250  °C  at  50  °C/min.  Thermal  desorber  parameters  were  a 
tube  desorption  temperature  of  250  °C  for  2.5  min  and  a  10  mL/min  split  flow.  The  tube 
purge/sweep  flow  was  50  mL/min,  the  trap  desorption  temperature  was  300  °C  for  2.0  min.  The 
sample  was  carried  into  a  split/splitless  injection  port  at  250  °C.  The  split  vent  was  turned  on  at 
0.5  min  with  a  purge  flow  of  50  mL/min  of  helium.  Retention  times  were  1  min  for  ethyl 
methylphosphonofluoridate,  2. 1  min  for  the  internal  standard  bromofluorobenzenc,  and  3.4  min 
for  3-hydroxymandelic  acid,  ethyl  ether  (impurity,  not  from  VX)  (Figure  1 ).  The  mass 
spectrometer  was  operated  in  Electron  Ionization  (El)  mode  and  scanned  from  35-300  amu  in 
2.78  s.  The  molecular  ion  for  ethyl  methylphosphonofluoridate  (m/z  =  126)  was  not  typically 
observed,  but  the  characteristic  fragment  ion  at  m/z  =  99  was  observed  (Figure  2). 

The  ability  of  the  Tenax  TA  tubes  to  collect  the  VX  vapor  was  \  crificd  by 
injecting  a  known  amount  of  VX  vapor  into  a  20  ft  copper  tube  that  had  flowing  air  and  showed 
that  the  VX  was  quantitatively  recovered. 
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Figure  1.  Retention  Times  of  Ethyl  Mcthylphosphonofluoridate  (MEPF),  1  min.  Internal 
Standard  Bromofluorobenzene  (BFB),  2.1  min.  and  Hydroxyniandelic  Acid.  Ethyl  Ether, 
3.4  min. 
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Figure  2.  Mass  Spectrum  of  the  I  1  min  Ethyl  Methylphosphonofluoridate  Peak. 


2.3  Analysis  of  Neat  Ton  Container  VX 

Two  samples  were  prepared  for  analyses  by  GC/MSD.  The  first  sample  was 
prepared  by  adding  2  pL  VX  to  200  pL  of  acetonitrile.  From  this  dilute  mixture,  a  10  pL  aliquot 
was  further  diluted  into  0.6  niL  of  methylene  chloride,  which  was  analyzed  by  GC/MS.  of  which 
1  pL  was  injected  onto  the  GC/MS.  The  compounds  of  interest  detected  were  2- 
(diisopropylamino)cthanol,  12.1  min;  diethyl  dimethylpyrophosphonatc,  16.2  min;  VX.  20.0  min; 
and  bis(diisopropylaminocthyl)  disulfide,  23.5  min. 

The  second  sample  was  prepared  by  adding  2  pL  BSTFA  and  2  pL  VX  to  200  pL 
acetonitrile.  The  mixture  was  then  heated  at  60  °C  for  20  min.  After  the  heating  period,  10  pL 
of  the  derivatized  mixture  were  added  to  0.6  niL  of  methylene  chloride  and  the  resulting  solution 
( 1  pL)  was  analyzed  by  GC/MS.  Compounds  of  interest  detected  were  trimethylsilyl  (TMS)- 
dcrivatized-ethyl  mcthylphosphonie  (EMPA)  acid.  1  1.5  min;  2-(diisopropylamino)cthanol. 

12.1  min;  TMS-derivatizcd-ethyl  mcthylphosphonothioic  (EMPT)  acid,  12.4  min; 
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S-tri methylsilyl-2-(diisopropylamino)ethanol,  16.0  min;  diethyl  dimethylpyrophosphonate, 

16.2  min;  VX,  20.0  min  and  bis(diisopropylaminoethyl)  disulfide  23.5  min. 

2.4  Experimental  Design  and  Data  Analysis 

The  data  were  analyzed  using  JMP  Statistical  Discovery  Software.  The  three 
variables  were  temperature,  drop  size,  and  air  flow  rate  at  three  levels  each.  Measuring  all 
combinations  of  these  levels  would  yield  27  conditions  (3  x  3  x  3);  the  cubic  composite  design 
chosen  required  nine  conditions,  which  can  be  described  as  the  vertices  of  a  cube  and  the  cube’s 
center.  This  collection  of  data  allowed  for  the  determination  of  the  major  contributing  variables 
and  interactions  among  variables,  although  in  this  particular  study,  not  all  of  the  data  points  in 
the  experimental  design  were  collected  due  to  the  length  of  time  of  each  experiment.  The 
substrate  temperature  (°C),  droplet  mass  (mg),  air  flow  (SLPM),  total  percent  VX  recovered,  and 
tunnel  identity  (four  similar  5  cm  tunnels,  named  a.  e,  k  and  1,  were  available)  were  treated  as 
variables  that  may  affect  the  raw  evaporation  rate.  Effect  interactions  among  droplet  mass,  air 
flow,  and  temperature  were  included  in  the  numerical  analysis. 


3.  RESULTS 

3.1  Time  to  Plateau 

The  VX  vapor  was  collected  until  the  concentration  of  agent  reached  a  plateau, 
often  ~1  x  10’* * * 4 5  mg/m  .  The  lowest  detected  vapor  concentration  was  2  x  10’6  mg/m3  (Figure  3). 
Only  the  50  °C,  0.92  pL,  405  SLPM  samples  consistently  gave  vapor  values  of  zero.  The  trend  in 
the  time  taken  to  reach  a  plateau  is  shown  in  the  cube  plot  in  Figure  4.  For  many  samples,  the 
concentration  at  this  time  was  near  the  short-term  exposure  limit  (STEL  required  by  the 
Occupational  Safety  and  Health  Administration  [OSH A])  of  1  x  10's  mg/m3.10  A  least  squares 
regression  (r*  =  0.94,  radf  =  0.91  )a  showed  that  the  significant  factors  were  temperature,  drop  size 
(in  mg),  air  flow,  %VX  recovered,  tunnel  ‘a’,  and  mass  *  air  flow  (Figure  5,  Table  1 ).  The 
equation  generated  was 

Time  to  plateau  =  6315  -  107*T  +  145  *  drop  mass  -  2.2  *  air  flow  -7.2  *  %VX  recovered  -  914 
*  tunncl[a]  -  0.42  *  j (drop  mass  -  4.3)  *  (air  flow  -  223)!  (1) 

3.2  Percentage  Recovery 

Seatterplots  of  the  variables  substrate  temperature  (°C),  droplet  mass  (mg),  and  air 
flow  (SLPM),  and  the  results  %VX  recovered  and  raw  evaporation  rate  showed  that  the  %VX 
recovered  and  raw  evaporation  rate  were  loosely  correlated  with  caeli  other  (r  =  0.64,b  Figure  6). 
None  of  the  other  parameters  were  correlated.  In  fact,  the  data  showed  that  the  %VX  recovered 
and  raw  evaporation  rate  data  were  distributed  throughout  their  respective  ranges  as  a  function  of 


J  f  estimates  the  proportion  of  the  variation  that  can  be  attributed  to  the  model  rather  than  to  random  error;  raJj‘ 

adjusts  for  models  that  have  different  numbers  of  parameters.  A  perfect  fil  is  r  =  1 .0.  JMP  Statistics  and  Graphics 

Guide,  Y'ersion  5,  by  SAS  Institute,  Cary,  NC.  2002.  pi 86. 

h  r  is  the  Pearson  Product-Moment  correlation;  a  perfect  fit  is  r  =  1 .0.  JMP  Statistics  and  Graphics  Guide,  Version 

5.  by  SAS  Institute.  Cary.  NC,  2002,  p376. 
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temperature,  drop  mass,  and  air  flow.  A  cube  plot  also  shows  the  ranges  in  %rceovcry 
(Figure  7).  The  percentage  of  agent  recovered  was  largely  random;  the  least  squares  regression 
had  a  low  correlation  coefficient  (r‘=  0.58,  ra(jfi=  0.37,  Figure  8). 

3.3  Evaporation  Rates 

The  evaporation  rate  was  calculated  by  summing  the  %VX  vapor  recovered, 
plotting  the  cumulative  %VX  loss  versus  time,  and  taking  the  slope  of  the  line,  which  was  the 
evaporation  rate  (Figure  9).  Examples  of  replicate  collections  of  evaporation  data  at  two 
different  conditions  demonstrate  the  degree  of  variability  observed  (Figure  10). 

Given  the  large  degree  of  variability  in  the  raw  evaporation  rate  (Table  2),  the  raw 
evaporation  rates  were  divided  by  the  %VX  recovered  to  yield  an  adjusted  evaporation  rate. 
Organizing  the  data  by  condition  and  taking  the  averages  showed  that  the  adjusted  evaporation 
rates  covered  a  much  smaller  range  than  the  raw  evaporation  rates,  and  major  trends  began  to 
emerge,  as  shown  in  a  cube  plot  (Figure  1 1,  Table  2). 

Numerical  analysis  of  the  data  was  performed  namely,  using  a  least  squares  fit 
for  the  raw  evaporation  rates  (mg/min)  as  a  function  of  substrate  temperature  (°C),  droplet  mass 
(mg),  air  flow  (SLPM),  %VX  recovered,  tunnel  identity  and  the  interacting  factors, 
temperature  *  mass,  temperature  *  air  flow  and  mass  *  air  flow(Figure  12,  Table  3).  The  r'  was 
0.90;  r;,dj  was  0.87,  and  the  statistically  significant  factors  were  temperature,  droplet  mass  (in 
mg),  %  agent  recovered,  temperature  *  air  flow  and  drop  mass  *  air  flow;  40  datapoints  were 
used.  Tunnels  ‘a’  and  ‘e’  were  significantly  different  from  tunnels  ‘k’  and  T;  and  an  equation 
was  derived  to  represent  the  relationship  (eq  2).  A  regression  that  included  the  %VX  recovered 
yielded  a  lower  r‘  of  0.89.  The  least  squares  equation  generated  was  then  used  to  predict  the 
evaporation  rate  for  each  sample  (Table  4). 

Raw  evaporation  rate  =  -8.08  x  10’4  +  1.41  x  10  4  *  T  +  4.05  x  10’4  *  Drop  mass 
+5.52  x  10'*  *  %VX  vapor  recovered 
+  1.30  x  10  *’*  (drop  mass-4.4)  *  (air  flow-207) 

+7.33  x  10’  *  (T-42.6)  *  (air  flow-207) 

+  1.08.x  10  '  *  tunnel[a]  -  9.83  x  1 0 4  *  tunnelfe]  (2) 

3.4  Combined  Effects 

The  simultaneous  effect  of  both  degradation  and  evaporation  was  calculated. 
Degradation  rates  for  VX  on  washed  glass  were  calculated  from  the  air-dried  sand  (ADS)  values 
at  22,  30,  40,  and  50  °C  and  the  washed  glass  (WG)  degradation  rates  at  22  °C  for  97%  or  purer 
VX.4  At  22  °C,  the  ratio  Awg/Aads  =  0.0059/0.0181  =  0.325.  The  rate  constants  at  50  °C  for 
ADS  (0.0912  hr’1)  and  the  interpolated  k  at  35  °C  (0.057  hr’1)  were  multiplied  by  0.325  to  gi\e 
an  estimate  for  the  degradation  rate  on  washed  glass.  Thus,  the  percentage  of  VX  that  had 
evaporated,  degraded,  and  remained  was  plotted  as  a  function  of  time  for  35  and  50  °C 
(Figure  13). 
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Figure  3.  Vapor  Concentrations  for  1  pL  Droplets  of  VX  Evaporating  from  Glass  at  35  °C, 
18  SLPM  (A, 111,0,0),  and  50  °C,  405  SLPM  (X, +,*-). 


2340 


900  - ► 

T  emperature 
35,  42,  50  °C 


Figure  4.  Average  Time  to  Reach  a  Plateau  in  the  Evaporation  Rate  of  VX  from  Glass. 
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Figure  5.  Actual  vs.  Predicted  Plots  for  the  Time  Taken  to  Reach  a  Plateau  for  the  Evaporation 
of  VX  from  Glass,  (a)  Least  Squares  Regression  Fit  (b)  Residuals. 


Table  1 .  Parameter  Estimates  for  Major  Effects  for  the  Time  to  Plateau  during 
the  Evaporation  of  VX  from  Glass. c 


Term 

Estimate 

Std  Error 

Prob> t 

Intercept 

63 1 5 

452 

<.0001 

Temperature,  T/°C 

-107 

1 1 

<.0001 

Drop  mass  /mg 

145 

31 

0.0001 

Air  Flow  rate/S l.PM 

-2.2 

0.7 

0.0056 

%VX  recovered 

-7.2 

2.9 

0.0229 

tunnelfa] 

-914 

285 

0.0044 

tunnelfel 

324 

198 

0.1 18 

tunnclfk] 

215 

141 

0.144 

(T-43.8)*(mg-4.3) 

-6.4 

4.4 

0.156 

(T-43.8)*(SLPM-223) 

-0.02 

0.10 

0.833 

(mg-4.3)*(SLPM-223) 

-0.42 

0.12 

0.0025 

c  Values  of  Prob>|l|  less  than  0.05  are  considered  significant. 
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Tabic  2.  Conditions  and  Experimental  Evaporation  Rates  for  the  Evaporation  of  VX  from  Glass. 


VX 

purity/ 

% 

n/a 

CD 

'c 

n/a 

© 

O' 

91.4 

DC 

OC 

n/a 

91.4 

n/a 

95.2 

95.2 

95.2 

n/a 

95.2 

90.4 

C3 

C 

n/a 

n/a 

n/a 

n/a 

n/a 

C3 

> 

DC 

c6 

Os 

93.8 

Average 

plateau 

time/ 

min 

1383 

2433 

344 

006 

3033 

Standard 
Deviation  raw 
evaporation 
rate/jig  min' 

© 

r- 

o 

© 

© 

0.2 

0.2 

Average  raw 
evaporation 
rate/pg  min  1 

0.6 

rj 

O 

2.7 

0.2 

0.9 

Standard 

Deviation 

%VX 

recovered 

22 

28 

r- 

Average 

%VX 

recovered 

sG 

22 

99 

r- 

33 

Time  to 
reach  a 
Plateau/ 
min 

1750 

O 

O 

oc 

1600 

2400 

ca 

2500 

2400 

na 

500 

420 

360 

370 

320 

250 

061 

1  100 

700 

3100 

na 

006c 

3100 

na 

na 

na 

Raw 

evaporation 
rat  c/ 
jig  min  1 

0.33 

0.71 

0.60 

0.15 

t  to 

0.12 

© 

o 

0.19 

4.30 

on 

3.20 

2.10 

2.10 

2.90 

ore 

r- 

o 

o 

0.36 

1.03 

0.75 

0.58 

1.10 

7.20 

3.70 

2.40 

%  VX 
recovered 

52.7 

45.1 

DC 

26.3 

24.1 

21.8 

37.1 

1 13.9 

33 

96.3 

53.8 

48 

58.3 

62 

r- 

27 

43.3 

13.3 

24.8 

49.3 

1  13.8 

61.8 

42.9 

<  -2  j 
u.  ^ 

r- 

DC 

c> 

oc 

181.3 

18.7 

18.7 

DC 

DC 

18.8 

r- 

DC 

405.6 

406.2 

405.4 

405.7 

405.6 

405.6 

405.6 

18.7 

18.7 

r- 

OC 

Os 

DC 

SO 

DC 

rn 

DC 

r- 

DC 

SO 

DC 

so 

DC 

Drop 

Mass/ 

mg 

0.92 

0.92 

0.92 

0.92 

0.92 

r  j 
O' 

o 

0.92 

0.92 

0.92 

C4 

O' 

© 

0.92 

0.92 

0.92 

n 

O' 

© 

rN 

O' 

o 

0.92 

0.92 

5.51 

5.51 

5.51 

5.51 

5.51 

5.51 

5.51 

Temp/ 

°C 

34.5 

34.9 

35.2 

35.3 

35.3 

35.5 

35.5 

35.8 

49.6 

49.8 

50.1 

50.1 

50.2 

50.3 

50.5 

50.4 

50.6 

34.5 

34.6 

34.8 

35.2 

41.4 

41.8 

41.8 

Code 

sC 

o 

A 

rn 

3a- 1 15 

3c- 1 57 

31-039 

r- 

O 

i 

O' 
r  t 

o 

A 

rn 

oeo-ie  1 

31-036 

CO 

(N 

i 

co 

uo 

O 

i 

CO 

31-033 

31-034 

3k-025 

3k-034 

3k-027 

^r 

(N 

O 

i 

c*D 

3k-026 

r\ 

n 

© 

A 

r^i 

3a- 1 16 

31-020 

DC 

6 

c^i 

3k-028 

3k-033 

31-032 
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Tabic  2.  Conditions  and  Experimental  Evaporation  Rates  for  the  Evaporation  of  VX  from  Glass  (Continued). 


VX 

purity/ 

% 

86.3 

86.3 

94.9 

88.3 

94.8 

80 

80 

89.2 

93.3 

93.3 

94.8 

n/a 

90.4 

n/a 

n/a 

90.4 

Average 

plateau 

time/ 

min 

1850 

2800 

2340 

717 

Standard 
Deviation  raw 
evaporation 
rate  pg  min  1 

ri 

- 

Cl 

rci 

Average  raw 
evaporation 
rate/pg  min 

C4 

Cl 

sC 

Standard 
Dev  iation 
%VX 
recovered 

27 

25 

24 

Average 

%VX 

recovered 

65 

57 

35 

43 

Time  to 
reach  a 
Plateau/ 
min 

1700 

2000 

na 

2400 

na 

3000 

3000 

o 

O' 
r  4 

2000 

2200 

3000 

1600 

§c 

700 

650 

Raw 

evaporation/ 
rate/ 
pg  min'1 

2.60 

3.20 

3.70 

3.90 

1.40 

2.50 

1.50 

0.39 

2.00 

2.70 

1.70 

0.26 

4.70 

5.00 

o 

p 

8.90 

°o  VX 
recovered 

58.6 

74.9 

40.3 

75.3 

53.5 

54.2 

43 

17.3 

42 

37.9 

29.9 

6.3 

79 

DC 

-r 

25.2 

roz. 

Air 

Flow/ 

SLPM 

181.6 

181.5 

181.7 

405.4 

405.4 

405 

404.1 

18.7 

18.7 

18.7 

DC 

DC 

18.7 

18.7 

405.4 

405.4 

405.5 

Drop 
Mass  / 
mg 

5.51 

5.51 

5.51 

8.27 

8.27 

r  j 

oc 

8.27 

r- 

Cl 

X 

r- 

ci 

dc 

r- 

Cs| 

DC 

8.27 

r- 

DC 

r- 

r4 

DC 

r- 

rt 

DC 

8.27 

ri 

DC 

c. 

£  o 

o  ° 
f- 

42.0 

42.2 

42.2 

34.7 

34.7 

34.9 

p 

cri 

rT- 

35.5 

50.7 

49.6 

50.1 

50.5 

50.6 

50. 1 

50.7 

50.6 

Code 

X 

|  31-037 

3  k -031 

31-027 

PC 

rE 

|  31-038 

|  3  k -03  9 

I  31-029 

cC 

i 

sC 

o 

1 

CT 

rn 

o 

1  31-022 

31-026 

|  31-021 

|  31-023 
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Figure  6.  Scatterplot  of  Temperature,  Drop  Mass,  Air  Flow  (SLPM),  %Recovery,  and  Raw 
Evaporation  Rate. 
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Figure  7.  Cube  Plot  of  %VX  Recovered  for  Evaporation  from  Glass. 
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Figure  8.  Least  Squares  Regression  Analysis  of%VX  Recovered  Data  as  a  Function  of 
Temperature,  Air  Flow,  Drop  Si/e,  Tunnel  Identity,  and  Time  to  Plateau. 


Figure  9.  Calculation  of  the  Evaporation  Rate  fora  6  pL  Droplet  of  VX  Evaporating  from  Glass 
at  35  °C  and  181  SLPM.  The  squares  show  the  data  that  were  used  to  generate  the  regression 
line  of  slope  0.001 1  mg/min. 
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Figure  10.  Cumulative  VX  Vapor  Recovered  for  1  pL  Droplets  of  VX  Evaporating  from  Glass 
at  35  °C,  18  SLPM  (A.D.O.O),  and  50  °C,  405  SLPM  ( X, +,*-).  The  solid  line  represents 
the  degradation  rate  of  VX  at  50  °C  and  the  dashed  line  represents  the  degradation  rate  of  VX  at 
35  °C. 
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Figure  1 1 .  Cube  Plot  of  the  Raw  Evaporation  Rates  (pg/min)  for  VX  Evaporating  from  Glass. 


20 


0  009- 

0.002- 

03 

— 

ro 

3 

O 

0  007- 

3 

■o 

0.0015- 

< 

QJ 

03 

cr 

— 

CO 

CD 

cr 

£ 

CD 

o.ooi- 

0.005-  .  ■ 

■ 

.  \ 

0.0005- 

Q_ 

03 

> 

0  003-  "  ■ 

_  j-  ■ 

cr 

o 

o- 

UJ 

CD 

£ 

03 

o.ooi-  r  ■  * 

> 

LU 

-0.0005- 

q: 

-  .  - 

> 

-0.001- 

-0.001- 

z> 

CTJ 

a: 

-0  002  .002  004  006  008  010 

-0  0015- 

Raw  Evap  Rate  Predicted  P<  0001  -0  002  .002  .004  006  .008  010 

RSq=0.90  RMSE=0.0007  Raw  Evap  Rate  Predicted 


Figure  12.  Actual  vs.  Predicted  Plots  for  the  Raw  Evaporation  Rates  for  VX  on  Glass,  (a)  Least 
squares  regression  tit  (b)  residuals. 


Table  3.  Parameter  Estimates  for  the  Major  Effects  Contributing  to  the  Raw  Evaporation  Rate 
for  the  Evaporation  of  VX  from  Glass/1 


Tenn 

Estimate 

Std  Error 

Prob>|t| 

Intercept 

-8.1  x  I03 

9.4  x  ur4 

<.0001 

Temperature,  T  /°C 

1.4  x  104 

2.3  x  10  s 

<.0001 

Drop  mass  /mg 

4.0  x  10'4 

5.3  x  1(T5 

<  0001 

Air  flow  /SLPM 

-1.0  x  10'6 

1 .0  x  1 0'6 

0.3004 

%VX  vapor  recovered 

5.5  x  10  - 

6.0  x  1 0'6 

<.0001 

(T-42.6)*(drop  mass-4.4) 

-8  x  10'7 

7.0  x  lO'6 

0.9106 

(drop  mass-4.4 )*(air  flow-207) 

1.3  x  10'6 

2.6  x  10' 7 

<.0001 

(T-42.6)*(air  flow-207) 

7.3  x  10'7 

1.5  x  10'7 

<.0001 

tunnelfa] 

1.1  x  10-’' 

4.5  x  lO"* 

0.0219 

tunnel[c] 

-9.8  x  lO”4 

4.3  x  10-4 

0.0316 

tunnelfk 

1.9  x  10’4 

2.6  x  10'4 

0.4716 

ll  Values  of  Prob>jt  less  than  0.05  are  considered  significant. 
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Figure  13.  Relative  Contributions  of  Degradation  and  Evaporation  for  the  Loss  of  VX  from 
Glass,  (a)  50  °C  (b)  35  °C. 
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Table  4.  Comparison  of  Predicted  and  Measured  Raw  Evaporation  Rates  for  the  Evaporation  of 
VX  from  Glass. 


Code 

Temperature/ 

°C 

Drop  mass/ 
mg 

Air  Flow/ 
SLPM 

Raw  Evaporation/ 
Rate 

Mg  min 

Predicted  Raw 
Evaporation  Rate/ 

Mg  min 

3k-016 

34.5 

0.92 

181.7 

0.33 

0.23 

3a- 1 15 

34.9 

0.92 

181.9 

0.71 

0.76 

3c- 157 

35.2 

0.92 

181.3 

0.60 

1.0 

31-039 

35.3 

0.92 

18.7 

0. 1 5 

0.22 

3k-037 

35.3 

0.92 

18.7 

0.1  1 

0.57 

3k-029 

35.5 

0.92 

18.8 

0.12 

0.45 

31-030 

35.5 

0.92 

18.8 

0.01 

-1.2 

31-036 

35.8 

0.92 

18.7 

0.19 

0.8 

3k-023 

49.6 

0.92 

405.6 

4.3 

5.3 

3k-035 

49.8 

0.92 

406.2 

1.1 

0.94 

31-033 

50.1 

0.92 

405.4 

3.2 

4.0 

31-034 

50. 1 

0.92 

405.7 

2.1 

1.7  1 

3k-025 

50.2 

0.92 

405.6 

2.1 

1.9 

3k-034 

50.3 

0.92 

405.6 

2.9 

2.5 

3k-027 

50.5 

0.92 

405.6 

3.1 

2.7 

3k-024 

50.4 

0.92 

18.7 

0.07 

-0.29 

3k-026 

50.6 

0.92 

18.7 

0.36 

0.81 

3k-022 

34.5 

5.51 

181.7 

1.03 

1.5 

3  a- 1 16 

34.6 

5.51 

1 8 1 .9 

0.75 

0.70 

31-020 

34.8 

5.51 

181.6 

0.58 

0.0 

3c- 1 58 

35.2 

5.51 

181.3 

1.1 

0.70 

3k-028 

41.4 

5.51 

181.7 

7.2 

6.2 

3k-033 

41.8 

5.51 

181.6 

3.7 

3.4 

31-032 

41.8 

5.51 

1 8 1 .6 

2.4 

1.8 

3k-038 

42 

5.51 

181.6 

2.6 

3.2 

31-037 

42.2 

5.51 

1 8 1 .5 

3.2 

3.7 

3k-03 1 

42.2 

5.51 

181.7 

3.7 

2.2 

31-027 

34.7 

8.27 

405.4 

3.9 

3.3 

31-031 

34.7 

8.27 

405.4 

1.4 

2.1 

31-038 

34.9 

8.27 

405 

2.5 

2.2 

3k-039 

35.6 

8.27 

404.1 

1.5 

2.3 

31-029 

35.5 

8.27 

18.7 

0.39 

1.0 

31-035 

50.7 

8.27 

18.7 

2.0 

2.3 

3k-036 

49.6 

8.27 

18.7 

2.7 

2.6 

3k-032 

50.1 

8.27 

18.8 

1.7 

2.1 

31-022 

50.5 

8.27 

18.7 

0.26 

0.37 

31-026 

50.6 

8.27 

18.7 

4.7 

4.4 

31-021 

50.1 

8.27 

405.4 

5.0 

5.4 

31-023 

50.7 

8.27 

405.4 

4.0 

5.1 

31-025 

50.6 

8.27 

405.5 

8.9 

7.5 
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4. 


DISCUSSION 


The  precise  reasons  for  the  large,  vexing  variation  in  %VX  recovered  and  raw 
evaporation  rate  remains  unknown.  The  raw  evaporation  rate  was  derived  from  the  cumulative 
%VX  recovered.  Thus,  the  loose  correlation  between  %VX  recovered  and  raw  evaporation  rate 
was  not  surprising.  The  least  squares  analysis  used  both  %VX  recovered  and  raw  evaporation 
rate  to  generate  the  predictive  empirical  equation  that  had  r  =  0.90. 

Since  the  drops  were  not  individually  weighed,  the  exact  mass  for  each  drop  was 
not  known,  which  would  also  add  variability  to  the  data.  In  addition,  small  changes  in  the 
amount  of  adventitious  water  present  in  the  VX  would  affect  the  VX  degradation  rate,  which 
would  in  turn  affect  the  %VX  recovered.  Independent  studies4  showed  that  the  VX  degraded 
with  a  half-life  of  1 6  hr  (960  min)  at  30  °C,  and  8  hr  (480  min)  at  40  and  50  °C.  Thus,  agent 
degradation  was  competing  with  agent  evaporation  in  these  experiments.  Comparisons  of  the 
measured  evaporation  rates  with  the  calculated  degradation  rates  indicated  that  both  processes 
occurred  on  a  similar  timescale.  For  the  50  °C  samples,  there  was  great  variation  in  the 
cumulative  %VX  values  (Figure  10),  although  evaporation  for  all  samples  had  ceased  at 
400  min,  indicating  that  all  the  VX  had  been  exhausted.  A  mid-range  sample  was  chosen  for  the 
comparison  of  the  two  processes.  The  calculated  comparisons  of  ev  aporation  and  degradation  at 
50  °C  in  the  bar  graph  indicated  that  VX  would  still  be  present  at  600  min,  but  the  vapor 
recovery  in  all  eases  indicated  that  the  evaporation  of  the  VX  had  ceased.  This  difference  may 
indicate  that  the  actual  degradation  rate  was  faster  than  the  rate  used  in  the  calculations,  possibly 
due  to  water  in  the  air  or  EMPA  in  the  VX.  In  addition,  as  the  VX  evaporated,  the  concentration 
of  EMPA  increased,  thus  increasing  the  degradation  rate  further  in  the  autoeatalytie  cycle.  As 
seen  in  Scheme  2.  EMPA  is  both  product  and  reactant  for  the  degradation  of  VX.  Thus, 
degradation  rates  in  an  open  system  may  be  faster  than  in  a  closed  system  at  the  same 
temperature.  At  35  °C,  the  calculations  indicated  that  VX  would  be  present  at  2400  min,  as 
observed.  The  degradation  rate  was  faster  than  the  evaporation  rate  at  35  °C,  whereas  the 
evaporation  rate  was  faster  at  50  °C.  The  rate  constants  for  degradation  did  not  consider  air  flow 
or  drop  size  effects,  and  thus  only  the  temperature  component  is  reflected  in  the  bar  graphs.  In 
other  words,  with  a  different  air  flow  or  drop  size,  the  evaporation/degradation  balance  will  shift. 

The  EMPA  formed  is  non-volatile;  thus  a  mass-only  method  (such  as  balances  or 
TGA)  for  following  the  evaporation  of  VX  would  be  difficult,  because  the  mass  loss  would  be 
due  to  both  VX  and  diisopropylaminocthyl  thiol  (DESH).  As  much  as  3%  EMPT  product  was 
detected  in  the  NMR  studies;  the  wind  tunnel  studies  were  not  designed  to  detect  the  expected 
concomitant  bis(2-diisopropylaminoethyl)  sulfide  product. 
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Degradation  Pathways  for  VX.’ 


5.  CONCLUSIONS 

During  the  first  1500  min,  the  VX  vapors  emanating  from  the  glass  substrate  were 
above  the  Immediately  Dangerous  to  Life  and  Health  (1DLH)  limit  of  3  x  10  '  mg/m’  (1DLH  is 
the  Occupational  Safety  and  Health  Admininstration  [OSHAj  limit  and  the  short-term  exposure 
limit  [STEL]  of  1  x  10’5  mg/m'). 10  Both  degradation  and  evaporation  were  equally  important 
contributors  to  the  concentration  decrease  of  VX  vapors  in  these  experiments.  The  exact  point  at 
which  no  further  agent  evaporated  depended  on  the  temperature,  drop  size,  humidity  of  the  air. 
and  initial  purity  of  the  VX. 
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